STATISTICAL STABILITY OF LORENZ ATTRACTORS 
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Abstract. We consider the robust family of Lorenz attractors. These attractors are 
chaotic in the sense that they are transitive and have sensitive dependence on the initial 
conditions. Moreover, they support SRB measures whose ergodic basins cover a full 
Lcbcsguc measure subset of points in the topological basin of attraction. Here we prove 
that the SRB measures depend continuously on the dynamics in the weak* topology. 
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1. Introduction 

The theory of Dynamical Systems initiated by Poincare's work on the three-body prob- 
lem of celestial mechanics and it studies processes which are evolving in time. The descrip- 
tion of the processes is given in terms of flows when the time is continuous or iterations of 
maps when the time is discrete. An orbit is a time-order collection of states of the system 
starting from a specific state applying the flow or the map. The main goals of this theory 
are to describe the typical behavior of orbits as time goes to infinity, and to understand 
how this behavior changes when we perturb the system or to which extent it is stable. In 
this work we are concerned with the stability of a system. 

Ergodic Theory deals with measure preserving processes in a measure space. One in 
particular tries to describe the average time spent by typical orbits in different regions of 
the phase space. According to Birkoff 's Ergodic Theorem, such times are well defined for 
almost all points, with respect to any invariant probability measure. However, the notion 
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of typical orbit is usually meant in the sense of volume (Lebesgue measure), which is not 
always an invariant measure. It is a fundamental open problem to understand under which 
conditions the behavior of typical (with respect to Lebesgue measure) orbits is well defined 
from the statistical point of view. This problem can be precisely formulated by means of 
Sinai-Ruelle-Bowen (SRB) measures which were introduced by Sinai for Anosov diffeo- 
morphisms |Si72j and later extended by Ruelle and Bowen for Axiom A diffeomorphisms 
and flows |BR75] . 

Definition 1 (SRB measures). Let /i be an invariant Borel probability measure for a 
flow (X*) t on the Borel sets of a manifold M. The hasin of /i is the set of points x G M 
such that 

lim — / (/?(X*(x)) dt = [ (f (i/i, for any continuous </) : M — )■ M. (1) 
T^+oo T Jo J 

The measure fi is an SRB measure if its basin has positive Lebesgue measure. 

The notions of basin and SRB measure can easily be extended to discrete time dynamical 
systems, simply by replacing the integral by a time series in ([1]). 

A fairly good description of the statistical behavior of orbits can be given by an SRB 
measure in the sense that, for a "big" (meaning "positive volume") set of points, the time 
averages of a physical observable (a continuous function on the manifold) of the system 
is accomplished simply by integrating the observable with respect to SRB measure (space 
average). 

In trying to capture the persistence of the statistical properties of a dynamical system, 
Alves and Viana |AV02j proposed a notion, called statistical stability, which expresses the 
continuous variation of SRB measures as a function of the dynamical system. This is a kind 
of stability in the sense that the outcome of evaluating continuous functions along orbits 
does not change much under small perturbations of the system. This is what we might 
be observed in computer experiments, where typically the picture obtained by plotting an 
orbit seems to be independent of the starting point and truncation errors. 

Next we in introduce the notion of statistical stability for vector fields. 

Definition 2 (Statistical stability). Assume we have a family X of vector fields endowed 
with a topology, admitting a common trapping region U on which each X & X has a unique 
SRB measure ^x- We say that X is statistically stable (in U) if the map X 3 X \ — )■ fix 
is continuous, where in the space of probability measures we consider the weak* topology. 

Our goal in this work is to prove the statistical stability of a family of vector fields 
associated to the Lorenz equations. 

1.1. Lorenz equations. Lorenz |Lo63] studied numerically a vector field X defined by 
the system of equations 

X = a{y — x), 
y = bx — y ~ xz, 
z = xy — cz, 

for the parameters a = 10, 6 = 28 and c = 8/3. The following properties are well known 
for this vector field: 
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(1) 



X has a singularity at origin which DX{0) has real eigenvalues 
0<-\3^ 2.6 <\i^ 11.83 < -As ^ 22.83; 



(2) there is an open set U, trapping region, such that X*(t7) C U for all t > 0. The 
maximal invariant set in f/, A = nf>oX*(t/), is an attractor and the origin is the 
only singularity contained in U; 

(3) the divergence of X is negative: 



By Liouville' Formula, the flow of X contracts volume. Thus, A has zero volume. 

Lorenz found with his experimental computations that the flow is sensitive with respect 
to the initial conditions near the attractor, i.e. even a small initial error lead to enormous 
differences in the outcome. It was a challenging problem to give a rigorous mathematical 
proof for this experimental evidence. Tucker |Tu99j gave a computer assisted proof that the 
original Lorenz system indeed corresponds to a sensitive robustly transitive non-hyperbolic 
attractor containing a singularity. Moreover, he proved that the Lorenz equations define a 
dynamical system with the behavior of the geometric model introduced by Guckenheimer 
and Williams [GW79j that we describe next. 

1.2. Geometric model. The geometric model of the Lorenz attractor is given by a vector 
filed Xq which is linear in a neighborhood of the origin containing the cube {(x, y, z) : \x\ < 
1, \y\ < 1, \z\ < 1}. The real eigenvalues Ai, A2 and A3 of DXo{0) with the eigenvectors 
along the coordinate axis satisfy 

< -A3 < Ai < -A2. 
Consider the square on the top of the cube given by 



Let r be the intersection of S with the local stable manifold of the singularity. The segment 
r divides E in to two parts 

S+ = {(a;,!/, 1) e E : X > 0} and E" = {(x, 1) G E : x < 0}, 

The return map from each of these parts to the planes x = ±1 is 

(x,?/, 1) I — > (sgn(x),?/|x|'', IxH 

where s = — A3/A1 and r = — A2/A1. The images of E^ by this map are curveline triangles 
without the vertexes (±1, 0, 0) and every line segment in J-" = {x = const flE} except F 
is mapped to a segment in {z = constflS'^}. The time r which takes for each (x, y, 1) G E\F 
to reach is given by 



Now we suppose that the flow takes the triangles back to the E in a smooth way as it is 
shown in Figure [Tal The time needed to pass from the triangles to E is relatively short with 
respect to the linearized region and it is a bounded function. In particular, the dynamics 
in the linearized region dominates all estimates of expansion and contraction. We assume 
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divX = dx / dx + dy /dy + dz /dz 



(a + 1 + c) < 0. 



S = Ux,y, 1) : -- < X,?/ < - 




r{x,y, 1) = -—log |x 



that each hne segment in {z = const fl S*^} is be taken back to another hne segment in 
{x = const n S}. Hence, the resulting Poincare map form E \ F into S has the form 

P{x,y) = {f{x),g{x,y)), (2) 

for some / : / \ {0} I and g : I \ {0} x I ^ I, where / = [-|, 

When the flow is it can be shown that the holonomies along the leaves in J-' are C^+" 
|Ma87| IPT93] . Moreover, by a strong dissipative condition on the equilibrium, i.e. 

r > s + 2, where r = — — ^ and s = — — ^, (3) 

Ai Ai 

it can be shown that J-" is then a foliation by an application of 2-normal hyperbolicity 
of a flow |HPS77j . So the holonomies along the stable leaves are maps. In particular, 
the one- dimensional quotient map is away from the singularity. 
The one- dimensional map / is described in Figure [lb] and satisfies: 

(1) / has a discontinuity at x = with side limits /(O^) = —1/2 and /(0~) = 1/2; 

(2) / is differentiable on J \ {0} and f'{x) > for all x G / \ {0}; 

(3) the limit of f'{x) is infinity clS X approaches , 

(4) f"{x) > for X G [-i, 0) and f"{x) < for x G (0, i]. 




(A) Geometric Lorenz flow ^^"^^"^ 

Figure 1 

The map g in equation ([2]) is defined in such a way that \dyg\ < < from which it 
follows the uniform contraction of the foliation J-", i.e. there exists a constant C > such 
that for any given leaf 7 of the foliation and .^1, .^2 £ 7 and n > 1, 

dist(P"(ei),P"(6)) < C"A"dist(ei,6). 

An orbit starting from a point on S will follow the linear vector field until one of the 
triangles and then it follows the flow coming back to S, and so on. For many proposes, 
the study of the flow can be reduced to the study of the two dimensional Poincare return 
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map P, and moreover the dynamics of this map can sometimes be reduced to the study of 
one dimensional map /, which is called the Lorenz map. 

1.3. Statement of results. A crucial fact about the geometric Lorenz attractor is that 
it is robust, i.e. vector fields sufficiently close in the topology to the original one con- 
structed as above also have strange attractors. Indeed, there exist an open neighborhood U 
in containing the geometric Lorenz attractor A and an open neighborhood U of Xq in 
topology such that for all vector fields X the maximal invariant set Ax = r\t>oX*{U) 
is a transitive set which is invariant under the fiow of X. This is a consequence of the per- 
sistence of an invariant contracting foliation on the cross section E for X G W |AP10t 
Theorem 3.10]. 

The smoothness property is also persistent for all nearby fiows, since the strong 
dissipative condition (|3]) is an open condition in the topology. For a C^-close vector 
field X to Xq, the leaves of J^x are close to those of J-" and it follows that fx is 
close to /. Thus, there exits c G [— |, |] which play for fx the same role of for /, and 
the properties of / in section 11.21 are still valid for fx on a subinterval [—b, b] for some 
< b < ^ close to |. 

Definition 3. We define the family X of Lorenz vector fields as a neighborhood of Xq 
with the following properties: 

(1) for each X & X, the maximal forward invariant set Ax inside U is an attractor 
containing a hyperbolic singularity; 

(2) for each X G A", S is a cross-section for the fiow with a return time tx and a 
Poincare map Px] 

(3) for each X & X, the map Px admits a uniformly contracting invariant foliation 
J-x on S with projection along the leaves of J-x onto / given by a map ttx; 

(4) for each X E X, the map fx on the quotient space / by the leaves in J^x is 
piecewise expanding with two branches; moreover, fx{x) > a/2 except at the unique 
discontinuity point cx and lim^_j.^± f'x{x) = +oo. 

(5) there is some constant C > such that for each X E X 

rx(0<-Clog|7rx(0-cx|. 

Observe that as the length of / is equal to one, then |7rx(0 ~cx| < 1 for all X G X. For a 
detailed exposition on the properties of geometric Lorenz fiows see e.g. [APlOt Section 2.3]; 
see also [AV12t Equation (9)] for the last property. The main goal of this work is to prove 
the following result. 

Theorem A. Lorenz vector fields are statistically stable. 

2. Preliminaries 

Consider the family X of Lorenz vector fields as in Definition [31 We assume that for 
each X E X the derivative fx is monotonic on each branch. On the other hand, l//x 
is bounded because /x > 1. Therefore l//x is monotonic and bounded and hence is of 
bounded variation. It follows from [Vi97a[ Corollary 3.4] that each fx admits a unique 
ergodic invariant probability fix which is absolutely continuous with respect to Lebesgue 
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measure A, whose density dfix/dX is a bounded variation function and, in particular, it is 
bounded. 

We point out that statistical stability results for piecewise expanding maps have been 
obtained in |Ke82] . According to |Ke82[ Corollary 14.] or |B(;97[ Theorem 11.2.2.], the 
family fx with X G A" is in the conditions of the results by Keller. Moreover, the density 
dJix/dX can be obtained by mean of a Lasota-Yorke inequality whose constants can be 
taken the same for all Lorenz maps; see [Vi97al Proposition 3.1.]. Therefore the density 
functions dfix/dX are uniformly bounded |Vi97at Corollary 3.4]. Hence we have: 

Proposition 2.1. Each fx for X & X is strongly statistical stable, i.e. fx ' — > dfix/dX is 
continuous with respect to L^-norm in the space of densities. Moreover, there exists M > 
such that for all X E X we have dfix /dX < M. 

By lifting up the measure fix, one may construct a unique P^-invariant probability 
measure fix on the cross-section E which is an SRB measure. For any bounded function 
: S -> M, we define 0* : J ^ M by 

0*(x)= inf 0(0, 

where ttx : S — )■ / is the canonical projection by stable leaves. The next result is proved 
in [APPVOQl Corollary 6.2]. 

Lemma 2.2. There exists a unique probability Px-invariant measure fix on S such that 

fix{<P)= lim / i<PoP^)*dfix, 
for every continuous function : S — )■ M. 

Actually fix is an SRB measure for Px- The uniform contraction of the stable leaves 
implies that the forward time averages of any pair ^i, ^2 of points on a same stable leaf for 
continuous function : S — )■ M are equal 

^ n—l _ n— 1 

hm -5^0(Pi(6)) = lim -5^0(Pi(6)). 

j=0 j=0 

Hence the inverse image of the basin of fix under ttx is contained in the basin of fix- This 
shows that the basin of fix contains an entire strip of positive Lebesgue measure, because 
the basin of fix is a subset of positive Lebesgue measure. 

On the other hand, since the density dfix/dX is a bounded function, we conclude that 
the return time is integrable with respect to fix- Then we can saturate this measure along 
the flow to obtain a unique SRB measure ^x for the flow whose support is the attractor 
Kx whose ergodic basin covers a full Lebesgue measure subset of points on the topological 
basin of attraction; see |APPV09| Section 7]. 

Proposition 2.3. The flow of each X E X has a unique SRB measure fix given for any 
continuous map (f : U M. as 

fix{^) = ^—- / ip{X{i,t))dtdfix{0- 
f^x{rx) J Jo 
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3. Statistical stability for the Poincare maps 



Here we prove the statistical stability of the Poincare maps on the cross-section E, i.e. 
the SRB measures fix depend continuously on the vector fields. 

Let {Xn)n>i be any sequence in X converging to X G A* in the topology. To shorten 
notations, we use subindex n instead of X„ for n > 1 and no subindex instead of X. 

Let {/) : E ^ M be an arbitrary continuous function. 

Lemma 3.1. Given m > 1 and e > 0, there is Uq = nQ{m, e) such that for all n > Uq 

i(0op™r-(0op™ri dx<e. 

Proof. Given m > 1, we can write J |(0 o P^)* ~ (0 ° P^y\dX as the sum 

/ i(0op-r-(0op-rirfA+ / i(0op-r-(0op-rirfA, 



(4) 



m—l 



where P„ = < ^ r„ o P^ > > and = N{m) is some large number. Now, by the last 

I i=0 J 

property in the definition of the geometric Lorenz fiows and the fact that the leaves of J-'n 
are nearly vertical lines, there is some constant Ci > such that 



m— 1 



^ |{x e / : -Clog - c„| > N}\ 

1=0 
m—l 



KBn) < 

1=0 
m—l 



_N_ 

e c 



i=0 



As is bounded, the first integral in (j4]) can be made arbitrarily small, provided A^ is big 
enough. 

We now estimate the second integral in (j4]). Considering 



m—l 



1=0 



1=0 



> 1 



we easily have that the second integral in (jlj) is bounded by 



m — l 

ToPi<N+l 

i=0 



oP'^)*|rfA+ / |(0oP„"^)* - ((/.oP'^)*|rfA. 



Observe that A(A„) is small for large n, because for large n a point belongs to An only if 
it belongs to some small neighborhood of the (finite) set of discontinuity lines of P*". As 
is bounded, we have that the second term in the last integral above is small for sufficiently 
large n. 
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It remains to control 



{T.T=o^roP^<N+l} 



1(0 o^r 



Observe that the points in {Xli^o^ roP*<A^+l| must necessarily be out of a neighbor- 
hood of the discontinuity lines of the map P™. If n is sufficiently large, then the same holds 
for P^. This means that the return time associated to these maps is uniformly bounded 
for large n. Then, just by the continuous variation of trajectories in finite periods of time, 
we can make |(0 o P™)* — {(f) o P"*)*| small for large n. □ 

Lemma 3.2. For any m > 1 we have 



lim 



o P™)*rf/i, 



Proof. Given m G N, then 



oP-)*rf/i„- / (0oP™)*rf^ 



< 



oP™)*d^„- / (0oP™)*d/x„ 



+ 



o p'^ydfi 



Using that the density of /2„ converges to the density of /2 in the L^-norm, by Proposi- 
tion I2.H and the fact that </> is bounded, we easily see that the second term in the sum 
above tends to zero when n goes to oo. So, we are left to prove that the first term converges 
to zero when n goes to infinity. In fact, using the uniform boundedness of the densities in 
Proposition 12 . 1 1 we obtain 

I dfLn 



O P-)*rf/i, 



'oP^rdfin 



<mJ |(0oP- 



oP 



d\ 

oPn*\dX 



dX 



which, by Lemma I3.H can be made arbitrarily small for n sufficiently large. 
Proposition 3.3. lim fln{4') = 



□ 



Proof. The compactness of S implies that is uniformly continuous and, therefore, given 
e > there exists 6 > such that 



10(6) - 0(6)1 < e, for every 6, 6 e S with dist(ei, 6) < S. 



(5) 



As we know, the rate of the contraction of the stable foliation on S is uniform for all vector 
fields in X. So, the first return maps are uniformly contractive. In particular, given 6 > 
there exists niQ > such that for all n we have 



diam(P™)(7) < 6, for all -y E J^n and m > uiq. 



(6) 



Take arbitrary numbers mi, m2 with m2 > mi > mo- Given x G /, let 7 be the leaf in 
containing x and 7m2-mi be the leaf in containing P^2~'"i(7). We have 

{(f) o = inf 0|p™2(^) = inf (plprn.^p^^-^.^^^y 

As G 7m2-mi, we also have 

(0opr)*(r^-"^Ka:)) = mf0ip^n(,„^_^) 

Then, since 7m2-mi contains P™^~™^(7), it follows from (|5]) and ([6]) that 

1(0 o prn^) - (0 ° PD^ir^-'^Ha^))! < 



Knowing that J ° -P™^)*'^/^" = y (0 ° P^^T ° "^^dfin, because /i„ is an /„ invariant 
probability measure, we obtain 

Consequently, the sequence I / (</> o P™)*(i/i„ j is uniformly Cauchy, because mo does 

\>' / m,n 

not depend on n. Hence 

lim lim / (0 o P^)* = lim lim / (0 o P^)* c//i„. 



and so 



lim /i„(0) = lim lim / (0 o P™)* dfin = lim lim / (0 o Pf )* d/i, 

I— >oo n—^oo m— >oo / m— >oo n— >oo 



and by Lemma 13.21 

lim /(0 o PJ^'ydfin = / (0 o P")*ci/i. 

n^oo J J 

Taking limit with m — 00 we complete the proof, by definition of jl. □ 

4. Statistical stability for the flows 

Now we prove Theorem 1X1 Let be a sequence in X converging to X G in the 
topology. Using again shortened subindex notation as in Section [3l we need to prove that 
fin ^ iiii the weak* topology. 

Let ip : U — M be a continuous function. We have 







By adding and subtracting the term 



—— / / ^iX{U))dtdfi{0, 

Un[rn> J Jo 



we have \Hn{'^) — /^(v')! bounded by the sum of two terms 



1 1 




r{0 

\v{X{U))\dtdm- (7) 
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and 




MO 



First we prove the following lemma. 




^{x{i,t))dtdm 



(8) 



Lemma 4.1. lim 



Proof. The term | J Tn d^in — / r is bounded by 



Tn d^r. 



T dfir, 



+ 



T djln 



T dpi 



W* ~ 



The second term converges to zero as n goes to oo, since — /i and the set of disconti- 
nuity points of r has /i-measure zero [Bi68t Theorem 5.1]. The first term is less bounded 
by 



M 



T\dX. 



Take a small rectangle A containing the singular curve of r and let n be large enough 
such that the singular curve of r„ is also contained in the rectangle. Since r and r„ are 
integrable, the integral l^n — T\d\ can be made as small as we want by taking the 
rectangle A sufficiently small. Outside the rectangle A, we have finite period of times and 
then the difference |r„ — r| can be made small by considering large enough n. 

Lemma 14.11 implies that ([7j) converges to zero as n goes to infinity, since 



□ 



1 




r{0 



\^{x{^,t))\dtdm < 



1 



1 



W\\oo^^{r) 



The next lemma implies that ([8]) converges to zero when n goes to infinity. 



Lemma 4.2. lim 

Proof. The difference 





^(0 



^{x{i,t))dtdm- 



ip{Xn{^,t))dtdjln{0 - J ip{X{^,t))dtdjl{0 



is bounded by M times 

\^{X4U))-¥>iX{U))\dtdX{0 





MO 



\^{X4U))\dtdX{0. (9) 



t{0 



The first term in 

"riO 




^ can be written as 

\^{x^{^,t))-ip{x{^,t))\dtdm 



A Jo 



IviX^iU)) - v{X{U))\dtdX{0, 



where A is a small rectangle containing the singular curve of r. The first term is bounded by 
2||</3||oo TdX and since r is integrable, it can be made small by taking A a small rectangle. 
The return time r is bounded on A'^ and, if n is sufficiently large, the same holds for r„. 
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In fact, the return times are uniformly bounded for large n. Therefore, by the continuous 
variation of trajectories in finite periods of time, we can make \ip{Xn{^,t)) — ip{X{^,t))\ 
small for large n, which implies that the second term converges to zero as n goes to infinity. 
The second term in (Q is bounded by 

llV'lloo J \^ ~ ^n\d\. 

which as we saw in the proof of Lemma [4.11 goes to zero when n goes to oo. □ 
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